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^ : ABSTRACT 

^ I We study the rotational distortions of the vacuum dipole magnetic field in 

^ I the context of geometrical models of the radio emission from pulsars. We find 

CN I that at low altitudes the rotation defiects the local direction of the magnetic 

field by at most an angle of the order of r^, where r^ = r/Ric, r is the radial 
' distance of the radio emission and -Ric is the light cylinder radius. To the lowest 

Q I (ie. second) order in r^, this distortion is symmetrical with respect to the plane 

vn ' 

O ' distortion which is asymmetrical with respect to the (f2, fl) plane is third order in 

^ ■ 
O 



Tn- These results confirm the common assumption that the rotational sweepback 
has negligible effect on the position angle (PA) curve. We show, however, that 
Qh| the infiuence of the sweepback on the outer boundary of the open field line 

O ' region (open volume) is a much larger effect, of the order of r^^. The open 

^ ' volume is shifted backwards with respect to the rotation direction by an angle 

rri ' 1/2 

' 5ov ~ 0.2sinarn where a is the dipole inclination with respect to the rotation 

axis. The associated phase shift of the pulse profile A0ov ~ ^.2rl!'^ can easily 
^ ' exceed the shift due to combined effects of aberration and propagation time 

■ delays (~ 2rn). This strongly affects the misalignment of the center of the PA 

curve and the center of the pulse profile, thereby modifying the delay-radius 
relation. Contrary to intuition, the effect of sweepback dominates over other 
effects when emission occurs at low altitudes. For rj, < 3 ■ 10^'^ the shift becomes 
negative, ie. the center of the position angle curve precedes the profile center. 
With the sweepback effect included, the modified delay-radius relation predicts 
larger emission radii and is in much better agreement with the other methods of 
determining Tq. 

Subject headings: pulsars: general — polarization — radiation mechanisms: non- 
thermal 
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1. Introduction 

There are independent observational arguments which imply that the pulsar radio emis- 
sion occurs in a form of a narrow beam centered (or roughly centered) on the magnetic 
dipole axis. In many cases the position angle (PA) of the observed linearly polarized ra- 
diation changes its direction by nearly 180° when our line of sight crosses the radio beam 
(eg. Lyne & Manchester 1988). If associated with the direction of B, this change of PA 
can be naturally interpreted as a result of our line of sight passing near the magnetic pole 
(Radhakrishnan & Cooke 1969). Moreover, the width p of the radio beam determined for 
different pulsars from the observed width of their radio pulse profiles scales with the rota- 
tion period P as the opening angle of the open field line region, ie. p oc P"^/^ (Rankin 1990; 
Rankin 1993). 

It is commonly believed that the emission region associated with the beam does not 
extend beyond the region of open field lines (hereafter called "open volume") which cross 
the light cylinder of radius Ric = c/Q (c is the speed of light and Q = Qz is the angular 
velocity of pulsar rotation). The angular size of the open volume at (small) radial distance 
r is equal to 9^^ ~ {r/ RicY^"^ and the cone formed by tangents to magnetic field lines at the 
rim of the open volume has angular radius of 9\, ~ 1.5^^^. The radial distance of the emission 
region has not been established so far: both a high-altitude emission region extending over 
a small fraction of ^q^, as well as a low- altitude emission region which fills in a much larger 
fraction of may be responsible for the same shape of the radio beam. 

The sweepback effect was first investigated in detail by Shitov (1983) who considered it 
to explain the observed dependence of radio luminosity of pulsars as a function of period. He 
estimated the magnitude of the rotational distortions of the magnetic field from the torque 
responsible for the observed slowing down of pulsars. He found that at moderate altitudes 
within the open volume, "near" the dipole axis, the direction of the distorted magnetic field 
defiects from the direction of the pure (ie. static shape) dipole barely by an angle 

5,b^l.2(^-^^ sin^a, (1) 

where a is the dipole inclination with respect to the rotation axis. Gil (1983) proposed 
the sweepback effect to explain why the separation between the main radio pulse and the 
interpulse observed in the profile of PSR B0950-I-08 is significantly different from 180°. 

In 1985 Shitov incorporated the sweepback effect into the model of pulsar position angle 
curves proposed by Radhakrishnan & Cooke (1969) and showed that the sweepback results 
in a lag of the profile center (measured as the midpoint between the outer edges of the 
pulse profile) with respect to the center, or the "infiection point" of the position angle curve. 
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Shitov emphasized that the lag of the profile center was a sum of two effects: not only the 
center of the PA curve is shifted toward earlier phases (with respect to the nondistortcd 
case) according to the cq. (1), but also the center of the open volume is displaced backwards, 
which contributes to the total effect. 

In most of subsequent investigations, however, the sweepback has been neglected, mainly 
on the basis of eq. (1). Blaskiewicz et al. (1991, hereafter BCW91) proposed a relativistic 
model of pulsar polarization which took into account two important effects overlooked by 
Shitov: the presence of the corotational acceleration and the aberration effect. An excellent 
result of their work was the "delay-radius" relation, according to which the center of the PA 
curve lags the profile center by 

A^scw ~ 4-^ rad, (2) 
-Kic 

where r is the radial distance of the radio emission. With no dependence on viewing geometry 
parameters (like the dipole inclination a, or the viewing angle (^hs between the rotation 
axis and the observer's line of sight), their relation appears to provide a powerful method of 
determining r. Equally important, the delay-radius relation depends neither on the observed 
width of the pulse profile W nor on the separation between the conal components in the 
pulse profile. Therefore, the altitudes of radio emission provided by eq. (2) may serve to 
determine which magnetic field lines are associated with the outer edge of the profiles and 
which field lines correspond to the maxima of conal components (Mitra & Rankin 2002; Dyks 
et al. 2004a). Von Hoensbroech & Xilouris (1997) used the delay-radius relation to probe 
the radius-to-frequency mapping at high radio frequencies. Given that the method is based 
on a measurement of tiny shifts (of magnitude usually being a small fraction of one degree) 
between the centers of the PA curve and the profile, it is extremely sensitive to the assumed 
geometry of the magnetic field. The latter was taken to be a dipole of static shape, with no 
rotational distortions. 

Gangadhara & Gupta (2001) proposed another relativistic method of estimating radio 
emission altitudes for pulsars with both core and conal components. By considering the 
effects of the aberration and the propagation time delays they showed that the core com- 
ponent lags in phase the midpoint between the maxima of conal components, if the core 
originates from lower altitudes than the cones, and if the cones are axially symmetric around 
the core in the reference frame corotating with the star (CF). Dyks et al. (2004a) revised 
their method and showed that the phase shift between the core component and the pairs of 
conal components is equal to 

A0BflH ~ 2-— rad (3) 

which provides another method for determining r without information about viewing geom- 
etry (nor W). Dyks et al. (2004a) used the above relation, and the results of work by Gupta 
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& Gangadhara (2003) to calculate r for 6 pulsars with well defined core-cone systems. As 
in the case of the delay-radius relation, the above formula holds only for the magnetic field 
which is symmetrical about the (Q, p) plane (where /2 is the magnetic moment of the pulsar 
magnetic field), at least as long as one associates the assumed symmetry of the core-cone 
system with the geometry of the underlying magnetic field. Any asymmetrical (with respect 
to the (r2, jl) plane) distortions of the magnetic field would pose a serious problems for the 
framework of the model leading to eq. (3). Again, based on eq. (1), any influence of the 
sweepback was neglected. 

Given that the above-mentioned methods of determining r arc so sensitive to the as- 
sumed symmetry of the magnetic field around Jl, it is important to have the symmetry 
hypothesis well justified. It is also important to revise this assumption in view of the unac- 
ceptably low values of r which are often being derived with the BCW91 method: as found 
in BCW91, the "delay radii" r^ei imphed by their method (eq. 2) are often smaller (in some 
cases by an order of magnitude — see fig. 29 in BCW91) than the geometrical radii r^eo 
determined with the traditional geometrical method based on the measurement of profile 
widths (Cordes 1978; Gil & Kijak 1993; Kijak & Gil 2002). This poses a real problem for the 
BCW91 method, because the geometrical radii, in the absence of strong refraction effects 
(Lyubarski & Petrova 1998), should be considered as lower limits of r (Dyks et al. 2004a). 
Although one could explain this disagreement in many different ways (eg. underestimated 
theoretical width of the open volume, systematically overestimated impact angles and dipole 
inclinations a etc.), we show below that the rotational distortions of the static shape dipole 
may account for a large part of the discrepancies between r^^i ^-nd Tgeo. 

Recently Kapoor & Shukre (2003) considered the aberration effect and the rotational 
sweepback to investigate the relative locations of core and cone components in the pulsar 
magnetosphere. Although included in the model, the sweepback is again estimated with the 
help of eq. (1). Being aware of the limitations of Shitov's estimate, the authors emphasized 
the need for derivation of a more advanced formula describing the rotational distortions of 
the magnetosphere. They noted that a proper derivation "should make use of at least the 
magnetic field given by the full Deutsch solution (Deutsch 1955)". 

Such an estimate based on the Deutsch solution was done by Arendt & Eilek (1998), 
who concluded that the rotation distorts the magnetic field by a magnitude of the order of 
r/R\c. Being much larger than the Shitov's estimate, this distortion would strongly affect 
results in BCW91, GG2001, Hibschman & Arons (2001, hereafter HA2001), and Dyks et 
al. (2004a). On the contrary, HA2001 noted that the leading terms in the difference between 
the Deutsch field and the rigidly rotating static-shape dipole are of the order of (r/i?ic)^. 
Recently, Mitra & Li (2004) emphasized that on the theoretical side there is a great need to 
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develop and understand the details of the sweepback effect. 

In this paper we investigate the rotational distortions of the pulsar magnetic field assum- 
ing the approximation of the vacuum magnetosphere. The twofold nature of the sweepback, 
first noticed by Shitov (1983) will be highhghted, and limitations in apphcability of eq. (1) 
will be clarified (Section 2). The significance of the sweepback for the relativistic model of 
pulsar polarization will appear to be much larger than previously thought, which will have 
serious consequencies for the delay-radius relation (eg. modification of eq. 2, Section 3). 



2. The rotational distortions of the dipolar magnetic field 

We follow previous investigators (Deutsch 1955; Shitov 1983; Barnard 1986; Romani & 
Yadigaroglu 1995, hereafter RY95; Cheng et al. 2000, hereafter CRZ2000) in assuming that 
the magnetic field surrounding the neutron star (NS) may be approximated by the vacuum 
rotating dipole. As in Barnard (1986), RY95, and CRZ2000, we assume that outside the 
NS the field is the same as of the star-ccntcrcd point dipole, ie. we neglect the near-surface 
modifications of the magnetic field by the conducting sphere of the neutron star, derived 
by Deutsch (1955) (cf. Yadigaroglu 1997). Hereafter, the magnetic field will be called a 
"retarded dipole" and will be denoted by ^ret- In Appendix A we give the cartesian and the 
spherical components of S^et (eqs. Al - Al and A15 - A17, respectively). 

We want to estimate how much the rotational sweepback distorts the magnetic field at 
low altitudes (r -C R\c) ■ One measure of this is the difference between the retarded magnetic 

— * — * — * 

field Sret and the magnetic field of the static-shape dipole B^t- The components of B^t can 
be calculated with the help of eqs. (A15 - A17) taken in the limit of Tn ^ 1 (ie. with the 
ratio Tn set equal to zero). We define the difference as: 

= 5ret - 5st. (4) 

In all formulae we assume that both the retarded dipole and the static-shape dipole are 
associated with the same magnetic moment /I, which at the time t = is in the (£, z) plane 
(time t is measured in the Lorentz frame in which the neutron star's center of mass is at 
rest). Thus, at any instant /2st = /?ret = A*, where 

/I = //(sin a cos x -\- sin a sin Q^t y -\- cos a z). (5) 

In the CF the components of i?ret at any point which corotates with the magnetosphere do not 
depend on t (cf. eq. A4). Therefore, one is allowed to choose any convenient value. We take 
i = (/I in the (x, z) plane) and constrain our discussion to the half of the magnetosphere 
with positive values of x. The positive values of the y coordinate then correspond to the 
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leading part of the magnetosphere and the negative y correspond to the traihng part. In 
cartesian coordinates and for i = the difference is: 



sm a 



1 



ABy — sin a 
ABz — sin a 



1 2 



-xzrl + Oirt) 



(6) 
(7) 
(8) 



In agreement with the remark of Hibschman & Arons (2001), the leading terms are second 
order in r^- The second order terms of AB^ and of ABz do not depend on y. The second 
order term of ABy is odd function of y {ABy{y) = —ABy{—y) + 0(r„)). These featm^es 
are important because the symmetry of any vector field B = {B^, By, B^) with respect to 
the plane requires the following relations to be satisfied: Bx{x,y,z) = Bx{x, —y, z), 



By{x,y,z) 



-By{x, 



-y,z), and Bz{x,y,z) = Bz{x, —y, z) (ie. the B^ and B^ components 



must be even, and the By component must be odd in y). 

The angle k between i?ret and i?st, to the order rj^, is given by: 



^^{firt + f2rl + fsrlY^' 



< 



(9) 



where the functions /i, /2, and /a, given in Appendix A, depend on x, y, z, and on the 
inclination angle a (but not on rn) and in general have magnitude of the order of 1, except 
from special locations in the magnetosphere which we discuss below. 

Thus, the rotation causes the magnetic field to deviate from Bgt by an angle which 
at most is second order in r^. Along the magnetic dipole axis^ of an orthogonal rotator, 
however, (ie. for a — 90° and (x, y, z) = (r, 0, 0)), one obtains fi — 0, f2 — and fs — 
so that K is third order in rn: 

'^^l^l (10) 

in partial agreement with the estimate of Shitov (1983). Beyond the orthogonal dipole axis, 
however, as well as at the dipole axis of non-orthogonal rotator, fiy^O and k may be much 
larger. On the {^,fl) plane /2 = (because /2 oc y, cf. eq. A5). Beyond the {fl,fl) plane 
(/i 7^ 0, /2 7^ 0), the first two terms in eq. (9) dominate and give: 



/2 



2/; 



1/2' 



(11) 



^By the "magnetic dipole axis" we understand the straight line containing the magnetic moment jl. In 
the case of the retarded dipole, the axis cannot be associated with any magnetic field line. 
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Let us estimate the angles ki and Kt for two points Pi{xi, yi, zi) and Pt{xt, yt, Zt) located 
symmetrically on both sides of the {Q, p) plane. Let us consider the particular case in which 
the points lie in the plane of rotational equator, close to the rim of the open volume of 
orthogonal rotator, ie. Xi = = r(l — r^^Y^"^, y\ = —y^ = rr}/'^, and z\ = = 0, with the 
positive value of the y coordinate corresponding to the point Pi on the leading side of the 
open volume, and the negative y for the trailing point Pt (cf. Fig. 1). Then eq. (A5) gives 
/i = 4-V„, /a = ±3-Vy^ which results in fa + S'V^ and Kt ^ - 

Both angles are considerably larger than the distortion of the magnetic axis given by eq. (10). 

Eqs. (6 - 8) imply, however, that up to the second order in rn, the rotational distortion 
of i?st is symmetrical with respect to the (fi, /I) plane (because Ai?^, is even, and the leading 
term of t^By is odd in y). This (approximate) symmetry, shown in Fig. 1, implies that 

— * — * 

beyond the dipole axis of the orthogonal rotator the angle k between B^-et and Bgt provides no 

— * — * 

information about the magnitude of the asymmetry of Bj-et with respect to the {Q, jl) plane. 

— * 

To estimate the asymmetry for points located beyond the (Q, jl) plane, one must there- 

— * 

fore use the difference between azimuths of B^^i at the points Pi and Pt: 



A0i_t = |0t| - 01 



Bret,y{Pt] 



Bret,x{Pt) 



^rct,y(Pi) _ sina 3 



— * — * 

where 0i is the azimuth of Pret at the point P\{x\,y\,Z]), 0t is the azimuth of Pret at 
Pt{xi, -yi, zi) (cf. Fig. 1), and 



R - ii 



X z ( X^ 

3 — cos a + I 3—5- — 1 1 sin a 

TV \ 



(13) 



Equation (12) clearly demonstrates that A0i_t is third order in rn, ie. the rotation induces 
the asymmetry of Pret with respect to the (f2, /I) plane with magnitude of the order of r\. 
Equation (12) is not useful in the immediate vicinity of the {^,11) plane (nor at the plane 
itself), because 0i changes sign to negative (ie. Pret is parallel to the {Q, Jl) plane) for locations 
with the tiny azimuth 

y _ 2 ^ f xz 1 2 3 



011 - - r„ h ^ ^ - K (14) 

— * 

(ie. within a narrow region on the leading side of the {^,jT) plane Pret.y < and eq. (12) 

— * 

gives a sum rather than a difference of azimuths). To estimate the asymmetry on the (Q, jl) 
plane one can use the difference of azimuths of Pgt and Pret at a given (the same) point: 

A0,_r ^ 0. - 0.et - 1^ - - /4r^ + |4|^ rl (15) 
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where the function /4(x, y, a) is given in Appendix A (eq. A14). Beyond the (Q, jT) plane 
the first term in this equation dominates (ie. A</)s_i. ^ fif^n) it is symmetrical with 
respect to the (f2,/I) plane (ie. odd in y, oc y). Therefore, just like k given by eq. (9), 
A(j)s-r provides no estimate of the rotational asymmetry there. On the {Q, pi) plane /4 = 
and A0s_r does measure the asymmetry which is of the order of r^. Along the dipole axis 
Bst,x = 2/isina/r^ so that A0s-r — 3~^r^, independent of a. 

Eqs. (9), (12), and (15) can be summarized as follows: the rotation changes the compo- 
nents and the direction of the dipolar magnetic field by ~ r'j^. To the order of this change, 
however, is symmetrical with respect to the (^2, p) plane. The asymmetrical change of the 
magnetic field direction has much smaller magnitude of the order of and is given by the 
second term of eq. (15). 

An immediate consequence of eq. (9) is that with accuracy of the order of rn, the rotation 
does not affect the shape of the position angle curve which depends on the local direction of 

— * 

Sret- This allows one to neglect the influence of the sweepback on the position angle curve 
as long as only the flrst order effects in r^ are considered. Given the tiny magnitude of the 
rotational asymmetry as defined by eq. (12), investigators often neglect its infiuence on the 
shape of the pulse profile as well (eg. BCW91, GG2001, HA2001). As we show below, this 
is not justified. 



2.1. The rotational distortions of the open field line region 

In the method of BCW91, the center of the pulse profile is most efficiently measured 
as a midpoint between the outer edges of the profile. Therefore, the method is based on 
the assumption that the outer boundary of the open volume is symmetrical with respect to 

— * 

the (17, p) plane. Due to the complexity of the magnetic field lines in the retarded case, we 
determine the outer boundary of the open volume numerically, by finding the magnetic field 
lines which are tangent to the light cylinder. The method is described in detail in Dyks et 
al. (2004b). 

Thick solid lines in Fig. 2 present the transverse shape of the open volume at low 
altitudes (rn ^ 1), calculated for the retarded magnetic field -Bret- More precisely, they 
represent the crossection of the outer boundary of the open volume with a sphere of radius 
r = 0.01i?ic centered at the neutron star. Different panels correspond to different dipole 
inclinations a. The magnetic moment p in all panels emerges perpendicularly from the page 
at the point (a;m,2/m) = (0,0). The thin circles have radius equal to r^^^ = rr^^, and are 
centered at the (0,0) point to guide the eye in assessing the asymmetry of the open volume 
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around p,. Given the small difference between the local direction of Bret and Bst (eq. 9), 

one may regard each panel to be permeated by the magnetic field of the static dipole with 
the straight magnetic field line emerging from the (0, 0) point toward the reader. The field 
is symmetric relative to the {Q, fl) plane, ie. with respect to the vertical line of = (to 
be imagined in each panel of Fig. 2). In the course of corotation, the contours of the open 
volume outer boundary move to the left in Fig. 2, ie. < correspond to the leading, and 
> to the trailing side. An observer's line of sight cuts the contours horizontally, moving 
left to right. ^ 

The following important conclusions can be drawn from Fig. 2: 1) the open volume is 
strongly asymmetric with respect to the {Q, pi) plane; 2) the magnitude of the asymmetry 
depends on y^a and thus on the impact angle /3; 3) regardless of the value of (3 the pulse 
window associated with the outer boundary of the open volume is always shifted backwards, 
ie. toward later phases; 4) the magnitude of this rotational asymmetry is very large, much 
larger than the local changes of the direction of the magnetic field caused by the rotation 
(eq. 9), and even larger than (for large a the "retarded contours" are (on average) shifted 

1 /2 

toward later phases by ~ 0.2^^^ — 0.2rn = 0.02 to be compared with rn — 0.01). 

This numerical result implies that none of the previous estimates (neither eqs. 12, 15, 
nor the Shitov's formula 1) provide a reliable measure of the rotational distortion of the 
open volume shape. The reason for this is that the boundary of the open volume is not only 
determined by the local (ie. low-altitude) direction of -Bret) but also (and most importantly) 
by the geometry of B^^t near the light cylinder, where ~ 1. At R\c all "higher order" effects 
become comparable in magnitude to the lowest order effects, in the sense that ~ 1 for any 
m. With the strength of the rotational distortions being very large at i?ic, a very different set 
of magnetic field lines is picked up as the last open field lines which form the boundary of the 
open volume. This "retarded" boundary is highly asymmetrical with respect to the (fl,//) 



^Note that our Fig. 2, calculated for f = in eqs. (Al - A3), differs significantly from analogous figures 

shown in Arendt & Eilck (1998) and in CRZ2000. The location of the retarded polar caps in their figures 
corresponds to the magnetic moment /x rotated by the angle flRns/c with respect to the (i, z) plane (ie. they 
assume t = Rns/c in eqs. Al -A3). At the same time, however, they assume t = {p. in {x,z) plane) to 
position the polar caps for the static case. Therefore, their figures do not inform us what is the relative 
position of the static and retarded caps in phase - only the caps' shapes can be compared. (To enable this, 
their retarded caps would have to be derotated by QR^s/c with respect to the static caps.) Also, note that 
the components of the magnetic field for the static-shape dipole given in CRZ2000 (eqs. A1-A3 therein) are 
for t = whereas the components for the retarded dipole (eqs. B2-B4 in CRZ2000) are for t = Rns/Ric- 
Their difference does not give eqs. (6-8): by overlooking this misalignment of the dipoles Arendt & Eilek 
(1998) incorrectly estimated that the rotational distortions of the magnetic field are of the order of Vn- For 
moderate dipole inclinations (a ~ 40° — 50°) the contours shown in Fig. 2 possess a notch rather than a 
discontinuous "glitch" suggested in Arendt & Eilek (1998) (see Dyks et al. 2004b for details). 
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plane (Fig. 2), and the magnitude of this asymmetry has httle to do with the low-altitude 
rotational distortions as estimated with eqs. (12), (15) and with Shitov's formula, because 

the low-altitude crossection of the open volume boundary is an image of the strong near-R\c 
distortions projected through the continuity of the magnetic field lines. Hereafter, we will 
refer to the asymmetrical distortion of the open volume with the terms "backward shift" or 
"displacement" of open volume. 

Another numerical result is that the contours shown in Fig. 2 are (with high accuracy) 
the same for any radial distance rn, as long as <^ 1, and as long as their size is normalized 
by r^^Qy = rrli"^ , as in Fig. 2. Contours calculated eg. for rn = 10^'^ or rn = 0.1 look exactly 
the same as those shown in Fig. 2. This means that for r^ the overall (ie. averaged over 
the impact angle /5) angular displacement of the open volume 5ov is a fixed fraction of the 
angular radius of the open volume ^^v- 

^ aei, « ar]!^ (16) 

with a ~ 0.2 sin a. The dependence a oc sin a has been determined by noting that 5o-v 
(measured from the star center) decreases with a (Fig. 2), whereas the corresponding phase 
shift A0OV 5oy/ sin a of the pulse profile does not (cf. Fig. 3, the discussion in Section 3.1, 
as well as eq. A8 in Dyks et al. 2004a). 

Equation (16) implies that the rotational displacement of the open volume has magni- 
tude comparable to the combined effects of aberration and propagation time delays (hereafter 
APT effects, of magnitude 2rn, cf. Dyks et al. 2004a) for r^ ~ a^/4 ~ 0.01, which is quite 
typical estimate of radio emission altitudes (eg. Gupta & Gangadhara 2003; Kijak & Gil 
2003). For rn ^ 0.01 the backward shift of the open volume dominates over the APT 
effects. For ^ 0.01 the APT effects dominate over the open volume shift. 



2.2. Twofold nature of the rotational distortions 

We find, therefore, that the nature of the rotational distortions at low altitude is twofold: 
in addition to the famous (but negligible) asymmetrical distortion of the local magnetic field 
direction (of magnitude ~ rf^, eqs. 12 and 15) the rotation shifts the open volume backward 

1/2 

by a much larger amount of ~ rn • The low-altitude shift of the open volume is not caused 
by the distortions of the shape of magnetic field lines at low altitudes — locally their shapes 
are pretty much the same as those of the static dipole. It is due to the strong near-i?ic 
distortions: magnetic field lines which are picked up as the last open magnetic field lines by 
the tangency condition at the light cylinder are located asymmetrically with respect to the 

— * 

(Q, jl) plane. The near symmetry of the field lines at low altitudes does not imply similar 
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symmetry of the open volume, because its shape is also governed by the geometry of the 
field near Ric. 

Although Shitov (1985) in his analysis of the phase shift between the pulse profile 
center and the center of the position angle curve neglected the APT effects, he did include 
the backward displacement of the open volume and did emphasize the twofold nature of the 
rotational distortions (cf. his fig. 2). However, he has not provided any simple estimate of the 
open volume displacement (like eq. 16). Therefore, in most studies following his work only 
the tiny local deflection of B has been considered (eqs. 12, 15 and Shitov's estimate), usually 
only to infer that the rotational distortions are negligible in generation of any asymmetry in 
a pulse profile. 



3. Implications for the relativistic model of pulsar polsirization 

Implications of the rotational displacement of the open volume for the relativistic model 

1 /2 

of pulsar polarization are profound, because this effect is lower order (~ rj , eq. 16) than 
the effects considered so far (aberration ~ rn, propagation time delays ~ rn). Although the 
PA curve is practically unaffected, the center of the pulse profile to which the PA refers is 
considerably displaced. 

Let us define a phase zero as a moment at which an observer detects a light signal emitted 
from the neutron star center when jl was in the (^1,/2) plane. As discussed in Dyks et 
al. 2004a, the total phase shift A$ between the center of the position angle curve and the 
pulse profile center can then be separated into two components: the shift of the center of the 
PA curve by A0pa ~ 2rn towards later phases with respect to the zero phase, and the shift 
of the pulse profile center toward earlier phases by A0pf with respect to the zero phase. Had 
the boundary of the open volume been symmetrical with respect to the {Q, jl) plane (as in 
the case of the static-shape dipole) the profile center would be shifted forward in phase by 
^0pf ~ — 2rn which would result in the total shift of Ar^ as initially predicted by BCW91. 
Due to the backward displacement of the open volume given by eq. (16), however, the forward 
profile shift is decreased by A0ov = Fr]/'^ (with F a/ sin a) so that A0pf = —{2rr, — Frn'^). 

Therefore, the delay-radius relation of BCW91 (eq. 2) becomes: 

A$ = A0PA + (-A0pf) « 2r„ + (2rn - Fr]/^) (17) 

with F in general being a complicated function of a, C and r^- The complicated form of 
F results from the complicated shape of the open volume boundary (Fig. 2) which implies 
nontrivial dependence of F on the impact angle and thereby on C,. Since the sign of the 
impact angle (5 — C, — a provides no information about whether the viewing trajectory is 
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poleward or equatorward (in the sense defined in Everett & Weisberg 2001) hereafter we will 
use the angle 

-{-(3, if « > 90° 

which is negative/positive for poleward/ equatorward viewing geometry regardless of whether 

1 /2 

a > 90° or not. Changes of F as a function of r = /9x/^b ~ Px/i^-^^n ) arc illustrated in 
Fig. 3 for the same angles a as those in Fig. 2. The functions F{t) were calculated for 
Tn = 0.01, however, they change little with rn, as long as rn <^ 1. Fig. 3 shows that F is 
confined to the rather limited range of 0.1 — 0.4 for any combinations of a and (. Also, F is 
always positive which implies that the displacement of the open volume results in a smaller 
phase shift A$ than predicted by the original delay-radius relation (eq. 2). The radio 
emission radii provided by the original delay-radius relation are therefore underestimated by 
a factor which may be very large for some parameters. 



3.1. The misalignment formula 

The delay-radius relation which includes the rotational distortions of the open volume 
(in the vacuum approximation) becomes: 

A$ fti 4rn - F{a, (, rn) r^^ (19) 

and it is plotted in Fig. 4 for a few values of F equal to 0.1, 0.2, 0.3, and 0.4 (solid, dotted, 
dashed, and dot-dashed curve, respectively). (For convenience, A$ on the vertical axis is 
in degrees, whereas cq. 19 gives A$ in radians.) The thick dashed line presents the original 
delay-radius relation of eq. 2. For rn < F^/16 (eg. for < 10~^ for F — 0.4, dot-dashed 
line), ie. for small emission radii, the phase shift A$ becomes negative (ie. the center of the 
PA curve precedes in phase the profile center), which is a new feature in comparison with the 
original delay-radius relation which always predicted the delay of the PA curve with respect 
to the profile. Since eq. (19) predicts that the center of the PA curve may either precede 
or lag the center of the profile, it will be referred to as the "misalignment" formula. For 
large radii rn ^ 10"^ the formula always predicts positive A$ and, for increasing rn, slowly 
converges to the original delay-radius relation. 

For any A$, the radio emission radii estimated with the original delay-radius relation 
underestimate those implied by eq. (19) (cf. Fig. 4). For A$ ~ 1°, the delay-radius relation 
underestimates r given by eq. (19) by a factor of 1.5 — 4 (depending on F). For A$ 0.3°, 
0.1°, and 0.01° the underestimate factor is in the range of 2 — 9, 3 — 22, and 16 — 220, 
respectively. For A$ fa the underestimation factor is formally infinite. In the absence 



-13- 



of effects described in Section 4, tfie negative values of the shift A$ cannot be lower than 
-FV16 > -0.57°. 

Unlike the original delay-radius relation, the misalignment formula predicts that the 
implied emission radii depend on the viewing geometry (especially on (5^ and that this 
dependence cannot be separated from the dependence on r-a'- even for fixed a and /3x, when 
the line of sight probes the magnetosphere deeper and deeper (ie. when rn decreases in 

1 /2 

Fig. 4), the absolute value of the parameter r = /9x/(l-5rn ) increases — the line of sight 
cuts through the open volume more peripherally. This departure from r = in Fig. 3 implies 
that the value of F changes (sometimes abruptly) with varying r-a, ie. for fixed a and the 
value of F is not fixed — it depends on 

The complicated behaviour of the misalignment formula (19), is exemplified in Fig. 5, 
which presents A$ as a function of calculated numerically for a = 45° and three values 
of Cobs = 43, 45, and 47° (circles, squares, and crosses, respectively). For large rn (and 
so for large A$) the line of sight crosses the open volume nearly centrally (/5x <^ 
T ~ 0), so that F ~ 0.1, regardless of the value of Cobs a-nd a (cf. Fig. 3). Therefore, all the 
three numerical solutions stay close to the analytical solution for F = 0.1 (thin solid line). 
For smaller rn, (and A$ < 1°), the numerical results diverge from each other: the case of 
Cobs = « = 45° (squares) remains close to the analytical solution with F ~ 0.1 (it follows 
eq. (19) with F ^ 1.07), because the parameter r oc /3x is fixed and equal to zero. 
In the case of Cobs = 47° (equatorward viewing, crosses in Fig. 5), the parameter r increases 
with decreasing because the line of sight traverses more peripherally through the open 
volume. As can be inferred from Fig. 3 (panel for a = 45°) this makes F increase through 0.2 
up to ~ 0.26 for T ~ 1, and accordingly, the numerical solution in Fig. 5 crosses the dotted 
line for F — 0.2 and approaches the vicinity of the dashed line for F — 0.3. At logrn ^ —3.33 
the line of sight just grazes the outer boundary of the open volume (r ~ 1.08). At smaller 
radial distances the line of sight does not penetrate the open volume. 

In the poleward case of Cobs = 43° (circles in Fig. 5), r becomes more negative with decreasing 
rn. Since the backward displacement of the open volume is stronger on its poleward side 
(cf. Figs. 2 and 3) the solution crosses the analytical curve for F = 0.2 (dotted line in Fig. 5) 
earlier than in the equatorward case (ie. at smaller |r|). At logrn — —2.43, the line of 
sight starts to cut the radiation beam above the notch visible in Fig. 2. This results in a 
discontinuous increase of F from ~ 0.25 up to ~ 0.39 (Fig. 3). Therefore, the numerical 
solution jumps to the vicinity of the dot-dashed line of F = 0.4 (Fig. 5). For more peripheral 
traverses (ie. for smaller rn and more negative r), F changes httle between 0.39 and 0.35, 
and the numerical solution departs only slightly from the F = 0.4 curve. For logrn ^ —3.25 
the line of sight misses the open volume. This minimum value of logrn differs slightly 
from the one for /3x — +2° (crosses) because in addition to the backward displacement. 
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the open volume is also slightly shifted (with respect to jl) towards the rotational equator 
(ie. downwards in Fig. 2). 



3.2. Determination of emission radius 

The above-described complicated form of the function A$(rn, a, Cobs) does not allow 
an easy determination of Tn. Numerical determination of Tn based on the known values of 
A$, a, and Cobs requires rather complicated calculations. Therefore, below we discuss the 
particular cases when easy derivation of rn is possible and then we propose a procedure, 
which allows the possible range of rn to be constrained in a general case of arbitrary viewing 
geometry. 

As can be inferred from Fig. 3, the value of F is (nearly) fixed and close to 0.1 whenever 

1 /2 

|t| <^ 0.4, ie. whenever \(3\ <^ 0.6rn . The condition is fulfilled for any r-^ and a when /? ~ 
(a ~ Cobs, squares in Fig. 5). The other case when the condition is fulfilled is when the 
measured shift between the position angle center and the profile center is large (> 1°) and 
the impact angle \(3\ is small: the radiation comes then from high-altitudes, where 6^^ ^ \(3\. 
This is why all numerical results shown in Fig. 5 approach the solid line of F = 0.1 when 
A$ increases above 1°. 

When F is fixed (like in the above-described cases), the equation (19) can be inverted 
to obtain analytical solutions for rn(A$): 



= 32 



-1 



8A$ + F^ + k (16A$F2 + F^) (20) 



where the parameter A; = I for A$ > whereas k = ±1 for A$ < 0. Thus, the solution for 
Tn is unique if A$ is positive, whereas for the negative A$ two solutions are possible. One 
can attempt to reject the smaller one of these solutions by using the theoretical constraints: 

rn > 1^ and r^>^ (3\ (21) 

The second constraint holds when the emission is limited to the open volume.^ In the two 
above-mentioned cases (when /3 or A$ > 1°), one can directly use eq. (20) with F — 0.1 
to calculate rn. 



■^In application to the retarded dipole field this condition is approximate, not exact. Therefore, the 
minimum values of logrn determined numerically for j3 = ±2° in Fig. 5 (circles and crosses) differ slightly 
from -3.266. 
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In a general case (including the common case when /3 7^ and the shift is small) one can 
constrain to a degree which depends on the information available. Without the knowledge 
of a and f3 one can obtain a rough estimate of based on the measured shift A$ by assuming 
some averaged, fixed value of F in eq. 19 (eg. F = 0.2). This is equivalent to igoring the 
dependence of F{a, r) shown in Fig. 3. For rn < 10^^ this may give results wrong by as 
much as an order of magnitude due to the large horizontal spread of as allowed by the 
range of F (Fig. 4). 

Therefore, a better approach is to use the following two-step procedure: 1) Using eq. (20) 
one can calculate the range of r^ allowed by 0.1 < F < 0.4 and by the conditions (21). Or, 
one can draw a horizontal line in Fig. 4 to determine graphically the range of r^ allowed by 
the range of F (the measurement error of A$ can easily be taken into account by drawing 
a horizontal strip instead of the line). For —0.0358° < A<l> < 0, two ranges of r^ are 
allowed and again the constraints (21) may be used to try to reject the lower range. For 
—0.573° < A$ < —0.0358°, the allowed range of r^ is between r""° and r™^^ where r™" is 
given by eq. (20) with k = —1 and F = 0.4 and r™^^ by eq. (20) with k = +1 and the same 
value of F = 0.4. 

2) With the knowledge of a and P one can more tightly constrain the allowed range of r^ 
by calculating Tmin = ^-^d Tmax = ''■('^JT™) where and r^^ were determined in 

the previous step. Then, from Fig. 3, using the appropriate panel, one can determine the 
narrower range of (Fmin, Fmax) corresponding to the range of (rmin, Tmax)- Then, one returns 
to the step 2) in which the tightened range (Fmin, Fmax) must be used instead of the original 
range of (0.1, 0.4). 

If the sign of (3 is unknown, the step 2) may (but does not have to) give two ranges of 
Tn. If the sign is positive (equatorward viewing) the range of F considered in the step 1) 
may be narrowed to 0.1 — 0.25 (Fig. 3). 

As an example we apply the method to PSR B0301+19 and B0525+21. In the case 
of B0301+19 (P = 1.38 s), BCW91 find A<I) = 0.2 ± 0.1° at 1.4 GHz. The range of r„ 
for F between 0.1 and 0.4 is {rf'Krlf-''') = (2.0 • 10-^ 1.2 ■ 10"^)^ if we allow for the 
error of A<I>, the range becomes (1.4 ■ 10"^, 1.25 ■ 10"^) (step 1). For = -0.96 ± 0.63° 
(Everett & Weisberg 2001, hereafter EW2001) Fig. 3 gives (rmin,rmax) = (-0.3, -0.1) (with 
the error of A$ included) or (rmin, Tmax) = (-0.5,-0.03) (including also the error of P^)- 
EW2001 derived a = 162.4 ± 11.8°, so that n - a ^ 17.6° and we may use the panel of 
Fig. 3 for q; = 20° to constrain the range of F to (F^m, -Fma^) = (0.14,0.22). This new 
range of F translates to (r™'', r™^) = (2.7 • 10"^, 4.6 • 10"=^) if the error of A$ is neglected 
(and the error of f3^ is included). Taking into account the error of A$ one finally obtains 
(r™°,r™^^) = (2.0 ■ 10^'^, 5.3 ■ 10~^). None of the constraints (21) narrows this range. This 
result is in agreement with the condition r^ > Tgeo/Ric- For the observed pulse width 
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W 15.9° (BCW91) and for a and (3 cited above, one finds r^^/Ri^ ^ 1.3 • 10"^ Our 
value of Tn = (2 — 5.3) • 10~^ implies that the emission associated with the profile edge 
must originate from magnetic field lines with magnetic colatitude 6*111 ~ (0.5 — 0.8)6'^^. The 
original delay-radius relation gives ^ 9 - 10~^ (3 times smaller) which results in 0^/6^^ > 1 
(emission from the region of closed field lines) . 

In the case of B0525+21 (P = 3.74 s), BCW91 find A$ = 0.3 ± 0.1° at 430 MHz. At 
this value of A$ the range of F = 0.1 - 0.4 translates into (rJJ"°, rJJ^^^) = (2.6 ■ 10'^, 1.2 ■ lO'^) 
or into (2 ■ 10"^ 1.3 ■ 10"^) jf ^^^g g^ror of A$ is included. For = 1.5 ± 0.08° (EW2001) we 
find (riiiiii^riiiax) = (0.15,0.39). EW2001 suggest n - a = 63.2° which implies that F ~ 0.1 
(Fig. 3, panel for a = 60°). Using this value in eq. (19) one finds that = (2.6 ± 0.6) ■ 10~^ 
which is two times larger than the value predicted by the original delay-radius relation 
(eq. 2). However, the value is still six times smaller than r^eo/Ric — 1-7 • 10~^ (calculated 
for W — 20.4°) and implies that the emission at the profile edge comes from magnetic field 
lines with colatitudes 9^ — (2.4 — 3.0)^^^, ie. from the closed line region. A large part of 
this discrepancy can be removed by considering different emission altitudes across the pulse 
profile. 



4. Altitude-dependent position angle swing 

When different parts of the pulse profile originate from different altitudes, the PA curve 
can no longer be described by the standard equation of Radhakrishnan & Cooke (1969). 
The simple analytical equation for the altitude-dependent PA swing is given in eq. (16) of 
BCW91. The formula predicts that if the radial distance r of the radio emission is uniform 
across the pulse profile, the entire PA curve is shifted rightwards (toward later phases) by 2r 
with respect to the zero phase. This effect is illustrated in Fig. 6: The thin solid line with 
dots is the PA curve for the fixed emission altitude of — 10~^, calculated numerically for 
the retarded dipole field with a — 45° and C, — 43°. In agreement with BCW91, its center 
lags the phase zero by 2rn rad = 1.14°, and there is no discernible sign of the sweepback 
effect (see Section 2). 

If the central parts of the pulse profile originate from much lower radial distance r than 
the edge, and if r <^ O.OlPic, the central parts of the PA curve do not exhibit any appreciable 
shift and nearly follow the undisturbed S-curve of Radhakrishnan & Cooke. The sohd fine 
which nearly passes through the (0, 0) point in Fig. 6b has been calculated for the emission 
from the last open magnetic field lines of the retarded dipole for the same a and C, as above. 
The corresponding radial distance of the radio emission as a function of phase is shown in 
panel a of Fig. 6. Within the central parts of the PA curve, the emission altitude is negligibly 
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small and the PA curve nearly follows the undisturbed S-swing (dotted line in Fig. 6b) . 

Thus, in the case when the central parts of the profile are emitted from much lower 
altitudes than the outer edge of the profile, the PA swing practically does not undergo the 
delay by 2rn. The misalignment formula (19) then becomes 

A$ fti 2r„ - Frl/^ (22) 

and is illustrated in Fig. 7 with the same layout as Fig. 4. One can see that in the considered 
case, the shift is much smaller than the shift predicted by the standard delay-radius formula 
(eq. 2, thick dashed line in Fig. 7). For a typical phase shift A$ ~ 0.3° eq. (2) underestimates 
r by an order of magnitude. The range of altitudes for which the shift is negative is much 
larger than in the case with the constant emission altitude (eq. 19, Fig. 4). 

In the considered case of the low-altitude emission within the profile center, the deter- 
mination of the radial distance of radio emission r-^ is performed as before (Section 3.2) with 
the only difference in that eq. (20) should be replaced by the inverse of eq. (22), ie. Fig. 7 
should be used instead of Fig. 4. Even allowing for the low-altitude origin of the central parts 
of the pulse profile, for B0525+21 at 430 MHz we still have < Vgco/Ric- Apparently, either 
the radio beam of this pulsar does not fill in the open volume or other efi'ects are important. 
These may include the overestimate of dipole inclination and/or the broadening of the 
observed pulse width due to the low energy of radio emitting electrons (all of which lead 
to the overestimate of rgeo) or refraction, eg. Lyubarski & Petrova (1998). However, our 
radial distance rn — 7 • 10~^ underestimates r^eol R\c by a factor of 2.5 whereas the original 
delay-radius formula underestimates it by a much larger factor of 13. 



5. Conclusions 

The rotational distortions of the vaccum dipole have twofold effect: in addition to the 
small changes of the local direction of the magnetic field, the region of the open magnetic 
field fines undergoes a strong distortion. 

— * 

The change of the local direction of B is second order in rn. We find, however, that 

— * — * 

it is symmetrical with respect to the (Q, jl) plane. The largest asymmetrical change of B 
direction is much smaller - third order in rn, in agreement with Shitov's estimate. 

The displacement of the open volume shifts the pulse window toward later phases with 
respect to the center of the position angle curve. The shift has the magnitude of the order 
ofr^l The open volume shift modifies the delay-radius relation if the center of the pulse 
profile is determined as the midpoint between the outer edges of the pulse profile, which are 



-18- 



assumed to lie close to the outer boundary of the open volume. At low altitudes, where effects 
of aberration and propagation time delays are small, the open volume shift dominates and 
may result in the center of the PA curve preceding the center of the profile (negative A$). 
A majority of pulsars exhibits positive A<l> which means that the radio emission altitudes 
typically exceed ~ 10~^i?ic. 

The radii derived with the misalignment formula exceed those derived with the original 
delay-radius relation by a factor which increases quickly with decreasing altitude (and A$). 
This explains the trend of the delay-radius relation to predict emission radii smaller than the 
geometrical radii. The underestimate may also be produced/enhanced by the low-altitude 
emission within the central parts of the pulse profile. When both these effects work together, 
the standard delay-radius relation may underestimate r by an order of magnitude even for 
the relatively large phase shifts (A$ ~ 0.3°). 

The influence of the open volume shift on the method based on the core-cone shift 
(eq. 3) is difficult to assess, because the locations of the conal maxima do not need to follow 
the open boundary of the open volume. If they did, eq. (3) would have to be replaced with 
^4>DRH ~ 2rn — Fr]/'^. The influence of the open volume shift on the geometrical method 
is small because the rotation increases the transverse size of the open volume insigniflcantly 
(by a factor smaller than ~ 1.2, cf. Fig. 2). 

The discussed distortion of the open volume is generated by a high-order effect which 
was being neglected for years on the basis that its "local magnitude" at low altitudes is 
small (~ r^). The actual importance of this effect appears to be much larger than that. 
This suggests that other "high-order" effects, eg. the longitudinal polar cap currents, of 

3/2 

magnitude ~ r'n , or the inertia of electrons, may be a lot more significant than their low 
altitude magnitude suggests. The toroidal currents due to the corotation of the charge-filled 
magnctosphcre (~ r^^J have been shown to notably modify the shape of the open volume, 
however, in a way which is symmetrical with respect to the (I) plane (cf. fig. 4.11 in 
Beskin et al. 1993). 

We thank U. Dyks for the derivation of eq. (9). JD thanks B. Rudak for all the years of 
fruitful collaboration. This work was performed while JD held a National Research Council 
Research Associateship Award at NASA/GSFC. This work was also supported by the grant 
... (JD) and by the NASA Astrophysics Theory Program (AH). 
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A. Reteirded vs static dipole 



In the reference frame with z || ^ the "retarded" field B^^x. of the vacuum magnetic point 
dipole with the magnetic moment given by eq. (5) has the following cartesian components: 



B 



ret,i; 



ret,j/ 



jsxz cos a + sin a ^ [(3^;^ ~ '^^) + 3xyrn + (r^ — x^)rn] cos(Qt — rn) + 
\ixy — (3a;^ — r^)rn — xyr^ sin(f2t — rn) j | (Al) 
[^yzQ,o%a + sina ^ [3xy + (3y^ — r^)rn — xyr^ cos(Qt — rn) + 

[(3^2 _ ^2) _ 3^^^^ ^ ^^2 _ ^2^)^2j g.^^^^ - ^n) ) } (A2) 

|(32;^ — r^) cos a + sin a ^ [3x2; + 3y2;rn — X2;r„] cos(ni — rn) + 
[3|/2; — 3a;2;rn — yzr-^ sin(f2t — rn) j | (A3) 

It can be derived by any of the methods described in Yadigaroglu (1997), Arendt & Eilek 
(1998), or CRZ2000. Also, it can be obtained by taking the limit of 7?ns/-Ric ^ 1 in the 
solution of Deutsch (1955). For any position ro) time to and time interval At it holds that 



B 



ret,z 



S,et (to + At, i?,(QAt)fo) = i?,(QAt)4et (^0, ro) , 



(A4) 



where i?^ (a; )r represents the rotation of the vector r by the angle io around the z axis. Thus, 
the dependence on the time t only reflects the rigid rotation of the pattern of distorted 
magnetic field lines around the z axis with the angular velocity Qi. The magnetic field Sgt 
of the static-shape dipole associated with the same magnetic moment /2 (eq. 5) is given by 
the same set of equations (Al - A3) with rn = 0. 

Using eqs. (Al - A3) with t = it can be immediately shown that the difference 

— * — * — * 

AS = Bret — Bst is given by eqs. (6-8). Using (4) and (6 - 8) one can find that the angle 

— * — * 

K = arccos(5ret • -Bst/(-Bret-Bst)) is given by eq. (9) with /i, /a, and fs given by 



and 
with 



fi^hi- gl, /2 = /i2 - '2-9192 
h = 29I + hz-9l- 2yi(^3 + hi) 



91 - -^c. 



92 



-2ci 



b—^ — 1 ) sin q; + b—^ cos a 



xy y 

— ;r sm q; H — 5- cos a 



(A5) 
(A6) 

(A7) 
(A8) 
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53 

hi 



X 



3 — 7— 1 sin q; — 7— cos o; 



= - Cl sin a 3— + 1 



X 



2 xy 
-- Cl sm a — 

3 r r 

' x"^ 5 



Cl sm q; 

8 



9 



sm a 



1 + 3 {xr~^ sin a + zr"^ cos of 
The function /4 in eq. (15) is equal to 



(A9) 
(AlO) 
(All) 
(A12) 
(A13) 



f =- (l^Y sing y / 



1 //i\2 sina y ( z , • A 

' ' ' 3- coscK + 4— smaj 



(A14) 



where B^t^x is given in eq. (13). 



The cartesian components of the retarded magnetic field given by eqs. (Al - A3) can be 

— * — * 

rewritten into the following spherical components of Sret in the reference frame with z || Q: 



B 



ret,r 



2fi 



{cos a cos 9 + sin a sin 9 [rn sin A + cos A] } 



— sin a [(r^ — l) sin \ + r^ cos A] , 

{cos q; sin ^ + sin a cos 9 [— rn sin A + (r„ — l) cos A] } 



where 



A = Tn + — ^t. 

Again, the (spherical) components of Bst are given by eqs. (A15 - A18) with = 0. 
Using eqs. (A15 - A17) in eq. (4) it can be easily shown that 



(A15) 

(A16) 
(A17) 

(A18) 



ABr 

AB. 



— sm a sm 9 



— sm a 



cos 5 — - sin 5 + (r„) 



12 

- - sin (5 cos 5 + (r^) 

2 3 



ABe — sin CK cos 9 



1 2 

- cos 5 rjl sin 5 + O (r^) 

2 3 



(A19) 
(A20) 
(A21) 



where 5 — (f) — Q,t. 
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Fig. 1. — The influence of rotation on the local direction of the magnetic fleld. The vectors 

— * — * 

of Bret and Bst are shown for two points in the equatorial plane of the orthogonal rotator 
(a = 90°). The points are located symmetrically with respect to the plane (the plane 

is orthogonal to the page and contains x). The rotation is to the right and it is assumed 
that X ^ y (region near the dipole axis) and x ^ Ric. With accuracy of r^, the retarded 
field Bret is symmetrical with respect to the (fi, /I) plane, ie. ki = Kt and |0t| = 4>i- More 
precisely, k,\ = Kt + 0{r^) and 0t = —<Pi — 0{r^), ie. nt < k,\ and |0t| > (p\- 
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Fig. 2. — Crossections of the open volume with the star-centered sphere of radius r ^ Ri^ 

— # 1/2 

calculated for the retarded dipole field -Bret (thick solid). The circles have radius of rrn 
and are centered on the magnetic moment jl which protrudes perpendicularly from the page 
at the center of each panel ((a;m,|/m) = (0,0)). The magnetic field i?ret which permeates 
each of the panels is nearly the same as that of the static-shape dipole with the axis parallel 
to p, (ie. it protrudes from the (0,0) point too). Rotation is to the left. The backward 
displacement of the retarded contours with respect to the (0, 0) points (or circles) results in 
the shift of the center of the pulse profile toward later phases. The shape of the contours 
does not depend on r and P as long as r <S it!ic (-P = 1 s, and — 0.01 was assumed in the 
figure). Their size scales as rvn^. 
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Fig. 3. — The parameter F of eq. 17 as a function of the parameter r = /3x/ (l-Srn ) for the 
same dipole inchnations a as in Fig. 2. Though the results were obtained for — 0.01, the 
curves change httle with as long as Tn <S 1. 
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Fig. 4. — The phase shift between the center of the PA curve and the profile center as a 
function of the radial distance of the radio emission. The four curves show the case with the 
rotational distortion of the open volume included (eq. 19) and correspond to different values 
of the parameter F shown in the lower right corner of the upper panel. The thick dashed 
line presents the original delay-radius relation (eq. 2) which does not include the swccpback 
effect. In the upper panel the positive shift is shown (the PA curve lags the profile). The 
lower panel is for the negative shift (the PA curve precedes the pulse profile) . Note that the 
original delay-radius relation significantly underestimates Tn, especially for small A$ and 
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Fig. 5. — Three numerical results for a = 45° (and different viewing angles C) overplotted 
on the curves from Fig. 4. The circles are for the poleward viewing geometry with ( = 43° 
(the curve is broken into two parts, one with positive, and the other with negative A$), the 
crosses are for the equatorward viewing with ( = 47°, and the squares are for a = ^ = 45°. 
For more details see text. 
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Fig. 6. — Influence of variations of emission altitudes across the pulse proflle on the shape 
of the PA curve. Thin solid line with dots (panel b) presents the position angle curve for 
emission from the fixed radial distance of = 0.01. Its center lags the phase zero by 
2rn rad 1.14°. The thick solid line in panel b is for the radio emission from the last open 
field lines, ie. it corresponds to different radial distances shown in panel a. Note that the PA 
curve for the case of varying r does not exhibit any noticeable lag. The spread in r visible 
in panel a) near = — 7° corresponds to the notch in the open volume which appears for 
moderate dipole inclinations (cf. Fig. 2, a = 45°). The dotted line (nearly overlapping with 
the thick solid) is the curve of Radhakrishnan & Cooke (1969), undisturbed by the special 
relativistic effects. 
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Fig. 7. — Dependence of the shift between the center of the PA curve and the profile center 
on the radial distance of the radio emission in the case when the central parts of the pulse 
profile originate from much lower altitude than its edge (ie. for the case marked in Fig. 6 
with the thick solid line). The layout is the same as in Fig. 4. Note the increased divergence 
from the original delay-radius relation (thick dashed line). 



